The Klein first integrals in an equilibrium system with 
electromagnetic, weak, strong and gravitational 

interactions 



Jorge A. Rueda a < b , R. Ruffini a < b < c , S.-S. Xue a ' b 

a Dipartimento di Fisica and ICRA, Sapienza Universita' di Roma 
P.le Aldo Moro 5, 1-00185 Rome, Italy 
b ICRANet, P.zza della Repubblica 10, 1-65122 Pescara, Italy 
c ICRANet, University of Nice- Sophia Antipolis, 28 Av. de Valrose, 06103 Nice Cedex 2, 

France 



Abstract 

The isothermal Tolman condition and the constancy of the Klein potentials orig- 
inally expressed for the sole gravitational interaction in a single fluid are here 
generalized to the case of a three quantum fermion fluid duly taking into ac- 
count the strong, electromagnetic, weak and gravitational interactions. The set 
of constitutive equations including the Einstcin-Maxwcll-Thomas-Fcrmi equa- 
tions as well as the ones corresponding to the strong interaction description are 
here presented in the most general relativistic isothermal case. This treatment 
represents an essential step to correctly formulate a self-consistent relativistic 
field theoretical approach of neutron stars. 

Keywords: Neutron Star Matter; General Relativity; Thermodynamics of 
General Relativistic Systems. 



1. Introduction 

The unsolved problems of supernovae theories as well as the necessity of pro- 
cesses leading to electrodynamical phenomena during the gravitational collapse 
to a black hole [l| lead to the necessity of critically reexamining the current 
treatment of neutron stars. In a series of articles (see 0,0]), we have recently 
developed the first steps towards a new consistent treatment for the descrip- 
tion of neutron stars, well beyond the traditional Tolman- Oppcnhcimer-Volkoff 
equations. 

First we have generalized the treatment of compressed atoms of Feynman, 
Metropolis and Teller to the relativistic regimes (see Q for details). There, it 
has been enforced self-consistently in a relativistic Thomas-Fermi equation, the 
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condition of /3-equilibrium extending the works of V. S. Popov Q, Ya. B. Zel- 
dovich and V. S. Popov H, A. B. Migdal et al. 0,0, J- Ferreirinho et al. @ and 
R. Ruffini and L. Stella Is] for heavy nuclei. Then, through the using of scaling 
laws, following [l(], 11 , 12 [, this treatment was extrapolated to compressed nu- 



clear matter cores at nuclear and supranuclear densities. Such cores have stellar 
dimensions and mass numbers A ~ (mpi anc k/m„) 3 ~ 10 57 or M core ~ M Q . In 
addition, they fulfill global but not local charge neutrality having electric fields 
on the core surface, increasing for decreasing values of the electron Fermi en- 
ergy reaching values much larger than the critical value E c = mlc 3 /(eh), for 
Ef = 0. The assumption of constant distribution of protons at nuclear densities 
simulates, in such a treatment, the confinement due to the strong interactions 
in the case of nuclei and heavy nuclei and due to both the gravitational field and 
strong interactions in the case of nuclear matter cores of stellar sizes at nuclear 
and supranuclear densities. 

In a subsequent work Q , we have generalized the above approach explicitly 
including the effects of the gravitational field by considering the most simplified 
nontrivial but rigorous treatment of a general relativistic system of neutrons, 
protons and electrons in /3-equilibrium. It has been there proved that the tradi- 
tional treatment for the description of neutron stars adopting the condition of 
local charge neutrality is not consistent with the Einstein-Maxwell equations and 
with microphysical conditions of equilibrium within quantum statistics. The role 
of the constancy of the general relativistic Fermi energy of each particle species 
pioneered by O. Klein [13| has been there emphasized and, the full system of 
equilibrium equations consisting of the Einstein-Maxwell and general relativistic 
Thomas-Fermi equations has been formulated. The corresponding solution of 
such a system of equations has been there given in the simplest possible exam- 
ple of a configuration of neutrons, protons and electrons in /3-equilibrium with 
electromagnetic, weak and gravitational interactions. New electrodynamic and 
general relativistic properties of the equilibrium configurations have been there 
illustrated. 

The aim of this article is to make an essential new step: we further proceed 
to the description of a system of neutrons, protons and electrons fulfilling strong, 
electromagnetic, weak and gravitational interactions. The essential role of the 
Klein first integrals is evidenced and their theoretical formulation is presented 
in the Einstein-Maxwell background. For the sake of generality the treatment is 
performed in the most general case in which finite temperature effects are also 
taking into account. We adopt throughout the article natural units H = c = 1. 



2. The Constitutive General Relativistic Equations 

The densities in the core of a neutron star exceed the nuclear density p nuc ~ 
2.7 x 10 14 g/cm 3 and may reach densities of order ~ 10 17 g/cm 3 at the verge of 
the gravitational collapse of the neutron star to a black hole. There is therefore 
the need of a consistent relativistic theory for the description of the interac- 
tions between the matter constituents. In particular, approaches for the nuclear 
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interaction between nucleons based on phenomenologic al p otentials and non- 
relativistic many-body theory become inapplicable (see 14, 15j|). 



A self-consistent relativistic and well-tested model for the nuclear interac- 
tions is the Walecka model (see [l|l[l3] for details). This model share common 



features with the model adopted by Bowers et al. in [ljj, |15| ; in both of them the 
nucleons interact through a Yukawa coupling and the flat spacetime has been 
considered to construct the equation of state of nuclear matter. The technique 
of constructing the equation of state assuming flat spacetime has been generally 
used since, as pointed out in [l4, 15 1, as long as p < 10 49 g/cm 3 the gravitational 



contributions to interactions between particles are negligible. However, when 
we turn to neutron star configurations at nuclear and supranuclear densities, 
it has been shown in Q how the solution of the Einstein-Maxwell system of 
equations is mandatory. 

In the often called extended version of the Walecka model, the strong inter- 
action between nucleons is described by the exchange of three virtual mesons: a 
is an isoscalar meson field providing the attractive long-range part of the nuclear 
force; uj is a massive vector field that models the repulsive short range and; p is 
the massive isovector field that takes account surface as well as isospin effects 



of nuclei (see also [13, [19( ) 



The total Lagrangian density of the system is given by 

C = C g + £ f + +£ P + C 1 + C int , (1) 
where the Lagrangian densities for the free-fields arc 

£ 7 = -\f^\ (3) 

C* = ^V^W-t/^), (4) 

£ w = —^O^ + ^mlu^, (5) 

£ P = —n^-RT + \m 2 p p^, (6) 

where f2 M „ = d^uj^ - d u u^, TZ^ = d^p v - d v p^ = d^A v - dyA^ are the 
field strength tensors for the a/ 1 , p and fields respectively, V M stands for 
covariant derivative and R is the Ricci scalar. We adopt the Lorenz gauge for 
the fields A^, lj^, and p M . The self- interaction scalar field potential U{a) is a 
quartic-order polynom for a rcnormalizable theory (see e.g. po|). The specific 
functional form of /7(a) is not relevant for the scope of this article, thus we will 
not adopt any particular form of it hereafter. 

The Lagrangian density for the three fermion species is 

£f = & (*7"£V - lfc> (7) 

i=e,N 
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where i^n is the nucleon isospin doublet, ip e is the electronic singlet, m, states 
for the mass of each particle-specie and D p = d p + T p , being the Dirac spin 
connections that satisfy the commutation relation 

hM = d^-T^ la , (8) 

where denotes the Christoffel symbols. 

The interacting part of the Lagrangian density is, in the minimal coupling 
assumption, given by 

Ant = -fjacripNipN - Sw^J^ - g p p p J^ + eA^Jif^ - eA^J^ N , (9) 

where the conserved currents are 

= Tp N ri>N, (io) 

= i>NT^N, (11) 
J£ e = fa-f^e, (12) 

J% = ^(j-Y^-fi>N. (is) 

The coupling constants of the a, to and p-fields are g ai g u and g p , and e is 
the fundamental electric charge. The Dirac matrices 7 M and the isospin Pauli 



matrices satisfy the Dirac algebra in curved spacetime (see e.g. [21 





= 




= 2.9^, 


{7 M ,7,} 


= 25* , 







(14) 
(15) 
(16) 
(17) 



The Einstcin-Maxwcll-Dirac system of equations is then given by 



(18) 
(19) 
(20) 
(21) 
(22) 
(23) 
(24) 

where the scalar density n s = ipNipN, the nucleon effective mass rh^ = tojv + 
g a a, and 

V& = g u u^ + g P rp» + e (^p) A», (25) 





G fin + 8irGT pl , 


= o, 






= o, 




2 V TV 

\-m u u - g^J^ 
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+ m 2 p P U ~ 9 P Jp 


= o, 


V^W + 


d a U(a) + g s n s 


= o, 




Vft) - m N ] ip N 


= o, 




-eA»)- m e ] ip e 


= o, 
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is the effective four potential of nucleons. The energy-momentum tensor of 
free-fields and free-fermions T^ v of the system (|3j)-([6j) is 



where 



rpjjLh 

1 



a 

rjif_LV 



rp(J,V 
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-FHF° 



(26) 

(27) 
(28) 

(29) 

(30) 



are the contribution due to free-fields and Tj V is the contribution of free- 
fermions which we discuss below. 



3. The Thermodynamic Laws and the Field Equations in the Spher- 
ically Symmetric Case 

We first introduce the non-rotating spherically symmetric spacetime metric 



Mr) dt 2 _ e X{r) dr 2 _ ^2 _ ^ ^2 ^2 



(31) 



where the v(r) and A(r) are only functions of the radial coordinate r. 

For very large number of fermions, we can adopt the mean-field approxima- 
tion in which fermion-field operators are replaced by their expectation values 
(see e.g. [17| for details). 

We write the nucleon doublet and the electronic spinor as 4>i = ^i{k)e~ lk ^ x ^ 
in the phase-space. Suppose that neutrons, protons and electrons, and the 
corresponding antiparticles, are in thermodynamic equilibrium with a finite 
temperature T. The occupation fcrmion- number operators of the "fc" -state, 
Afi(k) = ipj(k)ipi(k) with i = e,p, n, are replaced by their Fermi-distributions 



ft(k) = (^t(k)^t(k)) 



(32) 



where ks is the Boltzmann constant, fii and e^(fc) = \/k 2 + riif denote the 
single-particle chemical potential and energy-spectrum (we recall that for elec- 
trons 777 e = m e ). The sign '+' correspond to particles and ' — ' to antiparticles. 
We do not consider "real" bosons to be present in the system; the only distribu- 
tion functions involved in the computation are due to fermions and antifermions 
and therefore phenomena as Bose-Einstein condensation does not occur within 
this theory (see e.g. [IH for details). 
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It is worth to recall that all the thermodynamic quantities, e.g. k, e, T 
are written here in the local frame which is related to the coordinate frame by 
the Lorentz "boost" 

Ai a) =(tt a ,Xa,e a ,*a), (33) 

where u a = e"/ 2 5°, \a = e A / 2 <5^, 9 Q = r<5 2 , and $ Q = r sin 05^, being 5p the 
usual Kronecker delta symbol. 

The number-density n, of the z-specie, taking into account the antiparticle 
contribution is, within the mean-field approximation, given by 

d 3 k[f+(k)-ff(k)]. (34) 



(2tt) 



The contribution of free-fermions and antifcrmions to the energy-momentum 
tensor can be then written in the perfect fluid form (see e.g. pjj ) 

= {£ + V)u^u v - TV", (35) 

where it M is the four- velocity of the fluid which satisfies u^u^ = 1, and the 
energy-density £ and the pressure V are given by 

s= Si ' v = E ( 36 ) 

i—n,p,e i=n,p,e 

being Ei and V% the single fermion-antifermion fluid contributions 

Si = jAs J ' d 3 ke t (k)[f+(k) + fr(k)], (37) 

v< - |p^/^[^(*) + /r(*)]- (38) 

The equation of state (|55|) - (|55)l satisfies the thermodynamic law 

£+V-TS = n iMi> (39) 

where S = S/V is the entropy per unit volume (entropy density) and /Xj = 
d£ / dn.i is the free-chemical potential of the z-specie. At zero-temperature T = 
0, ^ = y/(Kf) 2 +fhf and n, = (Kf) 3 /(3n 2 ), where Kf denotes the Fermi 
momentum of the z-specie. 

The scalar density n s , within the mean- field approximation, is given by the 
following expectation value 

n s = = 7A3 E / d3fc 5fy i/^w + /r(*)i- ( 4 o) 

^ ' i—n,p ^ ! 

In the static case, only the temporal components of the covariant currents 
survive, i.e. {ip{x)^' l il){x)) = 0. Thus, by taking the expectation values of 



G 



Eqs. (fTU|) - ([l"3")) . we obtain the non-vanishing components of the currents 

Jq H = n ch u = (n p - n e )u , 
Jo = n b u = (n n + n p )u , 
Jfi = n 3 u = (rip - n n )u , 



(41) 
(42) 
(43) 



where rib, n pi nn and n e are the baryon, proton, neutron and electron number 
densities which are functions only of the spatial coordinates, and uq = iJgoo = 

Making a variation of Eq. (|39|) and using Eqs. (f36|) (f38|) and (|40"]) , we obtain 
the generalized Gibbs-Duhem relation 



i— n,p,e 



dV = riidfii — g a n s d(T + SdT, 

i 

which can be rewritten as 

dV = riidfii — g a n s da + \ £ + V — rij/i. 

i—n,p,e \ i—n,p,e 



(44) 




(45) 



where we have used Eq. (|39[) to eliminate S, and we have used the relation 
between the scalar density and the fluid energy-density 



dE 

drfiN ' 



(46) 



which follows from Eqs. AM])-® and (|4Uj) . 

Therefore, the Einstein-Maxwell equations (fT8|) — (f22 ]) . within the mean-field 
approximation, become 



-A(r) 



-A(r) 
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1 d 2 v~ 


[dr 


~ dr) 
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1 f dv 



2 \ dr 



dX 
dr 
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dr~) 
dX\ 
"dr) 
dX 
dr 



= -e A eJ c V 



c A [d a U(a)+g s n s 



[9pJ p - 



m 2 oj] 



m2 pP] 



(47) 
(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
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where we have introduced the notation luq = cj, po = p, and A = V. The 

metric function A is related to the mass M(r) and the electric field E{r) = 
_ e -(.v+X)/2y> througn 

e -, W = ! _ + Gr ^ 2(r) = ! _ « + , (54 ) 

where we have introduced also the conserved charge Q(r) = r 2 E(r). 

An important equation, although not independent of the Einstein-Maxwell 
equations (I47[) (1531) . is given the energy- momentum conservation law 

V^T"" = -g u J^ v - gpJpl^ + eJfF»\ (55) 

from which we have 

dV (S + V) dv da du dp dV 

IF = — Tr ~ 9 ° ns Tr ~ 9 ^ ~ 9 " J "Tr eJ «** ' (56) 

where we have used the energy- momentum tensor given by Eq. (|26[) . 

4. Constancy of the Klein potentials and /3-equilibrium 

Introducing the nucleon doublet and the electronic spinor in the wave-form 
ipi = ij)i(k)e~ lk i' x in phase-space, the Dirac equations (|24j) become 

( 7A1 /Cf -mi)i>i(k) -0, (57) 

where 

/Cf = fc^ 1 - Vf, V? = -eA". (58) 

In the mean-field approximation, making the quadrature of Dirac operators in 
Eq. ([57]) and averaging over all states "fc", we obtain the generalized chemi- 
cal potentials or, for short Klein potentials for electrons E ei neutrons E n and 
protons E p 

E e = - cV = c" /2 M e - eV, (59) 

E P = y/gooPp + + g P P + eV = c u/2 p p + g^us + g p p + eV, (60) 
E n = \/gooPn + g u u - g P P = c' y/2 /i„ + g^uj - g p p, (61) 

where we have used Eqs. ([l4) -((I7 |l and Eqs. ([32]), j34]), (p5]) - (f38]> . In the zero- 
temperature case, they are generalized Fermi energies for electrons E e = Ef , 
neutrons E„ = E% and protons E p = . 

Using the equations of motion for the fields p, lu and a, and using the gener- 
alized Gibbs-Duhem relation (|45[) , the energy- momentum conservation equation 
(l56l) can be rewritten as 



i—n,p,e ^ 

+ g^ribduj + g p n 3 dp + en ch dV — 0. (62) 
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The isothermal Tolman condition |23J (see also [13j) demands the constancy of 
the gravitationally red-shifted temperature 

^ + ^dv = 0, or e v/2 T = constant. (63) 



Such a condition can be used into Eq. (|62|) to obtain 

nid(c l ,/2 >;) + g u n b dw + g p n 3 dp + en ch dV = 0. (64) 

i— n,p,e 

Moreover, using the expressions (|59p - ([6"0)) of the generalized chemical potentials, 
Eq. (|64|) can be rewritten as 

J2 ntdEi = 0, (65) 

i—n,p,e 

which leads for independent and non-zero particle number densities Hi =/= to 
the constancy of the Klein potentials (f59 |) - (j6"Tj) for each particle-species, i.e. 



E e = e v/2 fi e - eV = constant, (66) 
E p = e" /2 fi p + V P = constant, (67) 
E n — e v ' 2 jjL n + V n — constant, (68) 



where 



V p = g^uj + g p p + eV, (69) 
V„ = g^uj - g p p. (70) 

In the case of nuclear matter in /3-equilibrium (assuming not trapped neu- 
trinos), the values of the constant Klein potentials (f5S]) — (f^5|) are linked by the 
condition 

E n = E p + E e , (71) 
which can be rewritten explicitly in terms of the chemical potentials as 

p n = P P + Pe + 2g p pe~ u/2 . (72) 



5. Concluding Remarks 

We have presented the self-consistent equations of equilibrium at finite tem- 
peratures for a system of neutrons, protons and electrons in /3-equilibrium within 
the framework of general relativity including quantum statistics, electro-weak, 
and strong interactions. In the mean-field approximation, we obtained the gen- 
eralized particle chemical potentials from the Dirac equations for nucleons and 
electrons. 

From the Einstein-Maxwell equations, the thermodynamic laws and energy- 
momentum conservation, we obtain the constancy of the Klein potential of each 
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particle-specie and of the gravitationally red-shifted temperature throughout the 
configuration, i.e. the first Klein integrals and the Tolman isothermal condition 
respectively. In the non-interacting degenerate case, following a minimization 
energy procedure, it was demonstrated that the thermodynamic equilibrium 
condition of constancy of the generalized particle Fermi energy of all particle 
species holds (sec E. Olson and M. Bailyn j24j). Such a procedure can be 
straightforwardly applied to the present case, being the final result given by the 
equilibrium conditions and (|67l) . 

The precise values of such constants are linked, in the case of nuclear matter 
in /3-equilibrium, by Eq. (|7ip , and their full determination needs the inclusion of 
additional constraints to the system, e.g. global charge neutrality (see e.g. Q). 

The correct implementation of such generalized Thomas-Fermi equilibrium 
conditions needs the self-consistent solution of the global problem of equilibrium 
of the configuration following from the solution of the Einstein-Maxwell equa- 
tions (|37|). (|4*5|) . (|50"|) - (j54")) . the general relativistic thermodynamic equilibrium 
conditions ([55)1 . (f6"6")l and (|57|) . together with the constraints, e.g. /3-cquilibrium 
and global charge neutrality. 

Thus, the full system of Einstcin-Maxwcll-Thomas-Fermi equations can be 
rewritten in the form 



r dr I 



• 



-\(t) 



1 ldv\ 1 

~2 + ~j 9 = -OirG^ , 

r z r dr I r^ 



v" + 
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d 2 UJ 
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d 2 p 
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E e = 


B V ' 2 p 
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E p = 


e* /2 Mp + 


E n = 


e"/ 2 Mn + 
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r(v' + A') 



4 

dv 
dr 
dv 
dr 
dv 
dr 



dX 
dr 
dX 
dr 
dX 
dr 



-47ree l//2 e A (n p - n e ), 



e A [d a U{a) + g s n s ] 



T 



V p = constant, 



V n = constant, 
constant, 



[dp ^ 



m 2 uA 



(73) 
(74) 
(75) 
(76) 
(77) 

(78) 

(79) 
(80) 
(81) 
(82) 



where the constants E n , E p and E e are linked by Eq. (|71|) and V Pi n is given by 
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Eq. (|S9"|). In particular, in the degenerate case T = 0, Eq. (fT5"j) becomes 



where V = eV+E e and we have used Eq. (f?9")) into Eq. ([75)1 . This equation is the 
general relativistic extension of the relativistic Thomas-Fermi equation recently 
introduced in Q for the study of compressed atoms. In addition, Eq. (|83[1 has 
been recently used to obtain the globally neutral configurations in the simpler 
case of degenerate neutrons, protons and electrons in /3-equilibrium (see Q for 
details). 

This paper was completed in honor to Oskar Klein in the occasion of the 
"weekly Oskar Klein Centre seminar series" given by one of us (R. R) at Stock- 
holm University. We would like to thank the referee for the constructive and 
very interesting suggestions. J.A.R. thanks D. Pugliese and S. Cacciatori for 
helpful discussions on the subject. 
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